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We study particle creation in the presence of bulk viscosity of cosmic fluid in the early uni¬ 
verse within the framework of open thermodynamical systems. Since the first-order theory of non¬ 
equilibrium thermodynamics is non-causal and unstable, we try to solve the bulk viscosity equation 
of the cosmic fluid with particle creation through the full causal theory. By adopting an appropriate 
function for particle creation rate of “Creation of Cold Dark Matter” model, we obtain analytical 
solutions which do not suffer from the initial singularity and are in agreement with equivalent solu¬ 
tions of ACDM model. We constrain the free parameter of particle creation in our model based on 
recent Planck data. It is also found that the inflationary solution is driven by bulk viscosity with 
or without particle creation. 


1. INTRODUCTION 

In the past twenty years, particle creation in cosmo¬ 
logical models has drawn the attention of many cosmol- 
ogists. In some works, particle creation is considered 
in the context of the thermodynamics of open systems. 
Within a FRW background, a description for matter cre¬ 
ation has been proposed by Prigogine and his coworkers 
[I|. They used the generalized form of the first law of 
thermodynamics for open systems to describe the flow 
of energy from the gravitational field to the matter field, 
resulting in the creation of particles. This is interpreted 
as an additional negative pressure, which leads to a re¬ 
interpretation of the energy-momentum tensor. It means 
that during an irreversible process, the expanding space- 
time causes a growth of entropy by producing matter. 
The negative pressure of particle creation might play 
the role of dark energy which explains the accelerating 
universe. Lima and Alcaniz examined the implications 
of some FRW cosmological models with matter creation 
against the observations @,01 and later it was shown that 
these models are consistent with the observational data 

i- 

To study the dissipative processes during the evolution 
of the universe such as entropy production, there is a ro¬ 
bust model called “Dissipative Cosmology” which is re¬ 
lated to scalar dissipation by the space-time symmetries 
which may be analyzed via the relativistic theory of bulk 
viscosity 0| . Dissipative phenomena in a homogenous 
and isotropic universe are related to bulk viscous pres¬ 
sure. Assuming some dissipative effects of bulk viscosity 
take place at the cosmological scale, one can explain the 
accelerating universe Q- Indeed, bulk pressures could be 
the consequence of either of two types of processes. In the 
First type, different components of the cosmic substra- 
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turn are coupled due to their different internal equations 
of state. In this case, cooling of sub-fluids with the expan¬ 
sion of the universe will be different because the system 
tends to move away from equilibrium. Therefore, it gives 
rise to a bulk viscous pressure ii- In the second type, 
one can suppose a perfect fluid whose viscous property 
comes into the picture during the cosmological particle 
production [13, jlj . One of the sources of an effective 
non-vanishing bulk pressure of the cosmic fluid may be 
the non-conserving particle pro duction which has been 
studied in some papers [H [fl, [13, [l3 ■ 

Some epochs in cosmology, like inflation and reheating 
(through dissipative processes) are far from equilibrium 
dynamics. So it is necessary to use relativistic theory 
of non-equilibrium thermodynamics for them. There are 
two types of dissipative models based on relativistic non¬ 
equilibrium thermodynamics. In the first-order t ype the¬ 
ory developed by Eckart, Landau and Lifshitz |14l . d, 
authors consider first-order deviation from equilibrium. 
Consequently, the solutions are non-causal and unstable. 
In fact, the neglected second-order terms are responsible 
for this problem. The second-order type is the theory by 
Israel and Stewart [Sd which removes the above prob¬ 
lem by including the neglected terms d- Although the 
equations of dissipative non-equilibrium thermodynam¬ 
ics are complex, but under a condition called “isentropic 
particle production” which is characterized by constant 
entropy per particle, one can drive a simple relationship 
between the viscous pressure and the particle production 
rate di d • Under this condition, the cosmic substra¬ 
tum is not a conventional dissipative fluid but a perfect 
fluid with varying particle number. 

In the present paper, we will exploit “isentropic” condi¬ 
tion to overcome complexity of solving the Israel-Stewart 
second-order equation in the framework of open systems 
and in the presence of particle creation. This paper is 
organized as follows. In section [21 we first review vis¬ 
cous FRW universe in the presence of particle creation. 
Although the results of first-order bulk viscosity model 
with particle creation are interesting, but as it is non- 
causal and unstable, it is not a satisfactory theory [2l|. 
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Therefore in section [3] we try to solve the Israel-Stewart 
causal equation under the ’’isentropic” particle creation 
condition, exploiting a scenario of “Creation of Cold 
Dark Matter” (CCDM) [23 - 1^ . We obtain non-singular 
analytical solutions for the inflationary and radiation- 
dominated epochs which are consistent with the equiva¬ 
lent solutions of ACDM model. We show also that the in¬ 
flationary solution can be a consequence of viscous pres¬ 
sure of cosmic fluid. Finally in section |4l the conclusions 
of the present work are summarized. 


2. VISCOUS FRW UNIVERSE WITH PARTICLE 
CREATION 

In “Dissipative Cosmology” with particle creation, the 
Einstein’s equation becomes 

~ j ( 1 ) 

where is the most general form of the energy- 

momentum tensor of cosmic fluid called “effective energy- 
momentum tensor”. This tensor can be assumed as 

= {P + P + Pc + - {p+Pc + U)g^,y, (2) 

where p, p, 11 and are the energy density, perfect 
fluid pressure, bulk viscous pressure and particle-frame 
4-velocity, respectively and pc is defined as the “creation 
pressure” of particles out of the gravitational field [l|, [H, 


{p + p)dN 

=- 

n dV 

where N is the particle number, V is the volume and n = 
N/V is the particle number density. Also the equation 
of state of cosmic fluid is supposed as p = (7 — l)p in 
which 0 < 7 < 2 is the “adiabatic index”. Accordingly, 
the continuity equation in the presence of bulk viscosity 
pressure and particle creation pressure takes the form 

p + 3{p + p + Pc + Il)H = 0, (4) 


where the dot denotes derivative with respect to the cos¬ 
mic time. 

In the second order theory of the non-equilibrium ther¬ 
modynamics, the particle 4-flux and entropy 4-flux are 
defined by = nu^ and — {tIP/ 2^T)u^, 

in which is the particle-frame 4-velocity, s is the spe¬ 
cific entropy, r > 0 the relaxation coefficient for transient 
bulk viscous effects, ^ > 0 the bulk viscosity coefficient 
and T > 0 the temperature of the fluid. The transport 
equation for the bulk viscous pressure is [iE[i3 


n -k rn = -^0 - -tU 


0 



t' 


T 


( 5 ) 


where 0 = 3H is the scalar expansion. The functionality 
of 5, r and T are determined according to phenomenolog¬ 
ical approaches. For example, from the equilibrium ther¬ 
modynamical e quat ion of state taken from lattice QCD 
calculations (^. l26| T = pp^ and ^ = wp with 

7 = 1.318, p = 0.718, r = 0.213, wq = 0.5 - 1.5 GeV, 
Tc ~ O.lQGeV and w = 8.2149 x 10“®. The viscosity 
coefficient ^ determines the magnitude of viscous stress 
relative to expansion in the limit r —> 0. In this limit we 
have the first order theory [Il[Il,i.e. S^^ = sm. 

The bulk viscous FRW cosmologies are mostly based 
on the non-causal and unstable first-order thermodynam¬ 
ics. This is the result of the complexity of second-order 
equations. Supposing the expression in the bracket to be 
very small, Eq. m reduces to 

n -k rri = -3CH. (6) 

One way of ensuring that viscous signals do not exceed 
the speed of light in the truncated theory is to adopt the 
phenomenological model t = ^/p Q. Assuming p (x 9'^ 
and 11/0^ <C 1, one can find 

n = -^H. (7) 

Putting Eqs. © and © and the equation of state into 
Eq. (U) and also using p' = 2pH'/H, one can easily 
obtain Hubble parameter as a function of scale factor 
[2lj. Although, this result may be interesting but, since it 
admits dissipative signal with superluminal velocities (it 
is non-causal) [l^ and also as all its equilibrium states are 
unstable [l^ it is therefore not a satisfactory relativistic 
theory. In the next section we try to solve Israel-Stewart 
causal equation in the presence of particle creation. 


3. SECOND ORDER THEORY AND CCDM 


For adiabatic creation, the particle creation rate is 
given by 


F 


iV 

N 


n R 

= -+3^- 

n R 


and therefore Eq. can be rewritten as 


( 8 ) 


Pc = - 


ip + p) y 
3H 


(9) 


From the Gibbs equation together with Eqs. © and 
([ 8 |). the variation of the entropy of the cosmic fluid can 
be written as 


nTs = -Ue-r{p + p). (10) 

Although there is entropy production due to particle cre¬ 
ation, but as it was explained in first section, one can get 
“isentropic particle production” with s = 0 and F > 0. 







3 


Therefore, by substituting Einstein’s field equations into 
Eq. (fTUl) with s = 0 one obtains [1, [H, [2^ 


r _ 2 H 

e~ 


( 11 ) 


In the present paper, we combine a scenario of “cre¬ 
ation of cold dark matter” (CCDM) with bulk viscous 
model. In this CCDM model (22l - [^ which is introduced 
as an alternative to ACDM model, a new kind of accel¬ 
erating flat universe with no dark energy which is fully 
dominated by cold dark matter is studied. The creation 
rate of these particles is given by [ 2 ^ 

r = (12) 


where a is a positive constant and fldm is the dark matter 
density parameter. One can assume that the cold dark 
matter (CDM) particles are described by a real scalar 
field maybe as remnant of incomplete reheating of infla- 
ton scalar field after inflation. In this model, the positive¬ 
ness of the dark matter density at high redshifts and in 
the matter-dominated era requires that a < 0 ^,, where 
flm = ^dm + ^baryon = 1 is the total density parameter 
in the flat FRW universe. However, we should constraint 
the value of /3 = for our model independently. In 

the following, by using Eqs. m and (HU we try to solve 
Eq. dni). 

Here we adopt a standard power-law form for the tem¬ 
perature T and bulk viscosity coefficient ^ as [1, [13, 


where k = and I = 

The analytical solution of Eq. (IT51) for a general “m” 
is complicated but in the case of m = 1 / 2 , it is solvable 
easily and we obtain 

where /3, = /3 — 1 and ’’Ci” is a constant of integration. 
As H is large enough for the inflationary epoch, one can 
neglect the second term of the above equation and con¬ 
sequently it converts to Eq. (7) in the first order theory. 
Conversely, in late time that H becomes very small, the 
second term dominates. Therefore, the second term of 
Eqs. (fl^ shows the difference between first-order and 
second-order approaches for the bulk viscous pressure. 
Using Eqs. o, (0 and m we get 


H' - 


-fiR) 

^ R 


^{R)Rj3^{k 1 2 )_ 


H 


3R 


( 20 ) 


Here, we use a simple one-parameter function for “adia¬ 
batic index”, 7 ( 7 ?) which can describe a transition from 
an inflationary to a radiation phase [ 1 ^ 


7(i?) 


4 A^{R/R^f + (a/2)(i?/i?,)“ 
3 A^{R/R^Y+ {RIR^Y 


( 21 ) 


T = fipY (13) 

e = (14) 


where p,, and “m” are non-negative constants and r = 
(7 — l)/ 7 . Using the Friedmann equation p = and 
Eqs. (ITSll and (14), the terms in the bracket of Eq. 0 
can be written as a function of Hubble parameter 


T HE H T , , 

-=2r-, j = -=2(m-l)-, (15) 


On substituting Eq. (1151) and H = Hli into Eq. 0 we 
obtain 


n‘ -hH 


1 

tH 


3H 

2h 


1 5 1 

(1 -b r) — 
\ -r 



(16) 


where H denotes derivative of H with respect to Hubble 
parameter. By using Eqs. (fTTl) and (14), Eq. (00 can 
be rewritten as 


H -bH 
-18iJ2 

7 (r - 377) ■ 


2 3 1 1 

- 1 - (1 -b r) — 

77r7(r - 377) 7(r - 377) ^ ’h 


(17) 


Exploiting Eq. (fT^ and also the phenomenological model 
T = E^jp m, we obtain 


n‘ -b [kK 


-2m 


-bZ77-0n = -6 


77 


7(/3- 1 )’ 


(18) 


where A, is a constant, 7?* is a reference value of R and 
0 < a < 1 is a free parameter which has a relationship 
with the power of the cosmic time t during the infla¬ 
tionary epoch. Substituting Eq. (|2T|) into Eq. (l20l) and 
integrating, we obtain 


2Ci 

+ 02 , 




- —4/^7 


(77? 




7?3//3.7 = 
- —4/57 


7 j( 3 // 3 . 7 )-I -1 

( 22 ) 


where S = 2-\/3fo/(2 + and C 2 is a constant of 

integration. 

For 7? <C 7?*, 7 ( 7 ?) —>• 2a/3 and the Hubble parameter 
is given by 


77(7?) = 


201 


12a/3* + 9(5 


,C 2 ( —)®/^7?9/2/3.^ 
7?» 


- 5/3*a/3 


(23) 

It is noticed that with or without particle creation, the 
bulk viscous coefficient gives us an inflationary solution 
and a positive energy density for initial vacuum as 7? —>■ 0. 
In this way, neglecting the terms involving (5 one can ob¬ 
tain a non-singular de-Sitter solution as 7? oc exp(77fc7), 
where 77^ = (20i/9-\/35o)~^““^'^- la this picture in¬ 
flation appears when the viscous stress becomes large 
enough to create sufficient negative effective pressure. 
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Comparing with the Hubble parameter value in the scalar 
field inflation, we can find the values of the constants 
during inflation by using Hk = with 

the primordial scalar amplitude ln(10^°Hs) = 3.089 [s^, 
where M^i is the Planck mass and ”r” is the tensor-to- 
scalar ratio. 

For R which leads to 7 (i?) —>■ 4/3 we find a 

solution for H 


H{R) = 




25/3./3 


(24) 

On the other hand, the Hubble parameter for the ACDM 
model in the radiation era (rad) and in the absence of 
curvature and cosmological constant (which is certainly 
an acceptable approximation for relativistic high redshift 
universe) is 


H{R) ~ 3 X 10-^^cm-\nradhY^^ ■ (25) 

where Vt^ad ~ 4.4 X 10 ® [3l| and i?o is the value 
of scale factor at the present time. When in the Eq. 
(IMl) . the particle creation dominates over bulk viscous 
coefficient (for large R) we obtain power-law solutions 
for Hubble parameter as H <x. R'^^*. In order to con¬ 
strain the free parameter /3 of the model and have an 
acceptable power-law form for H in radiation era, we 
need 0 < /3 < 1, i.e. 0 < a < ^dm- Considering the 
recent Planck data for cosmological parameters we find 
a < fldm = 0.308±0.012 [s^. Therefore, comparing Eqs. 
(1^ and (1^ with each other for /3 —>• 1/2, i.e. a —>■ 0.154 
one obtains H ^ R~^ and consequently putting away the 
first negligible term in the bracket of Eq. (l24l) we find 
C~^^^RllA^ = 1.99 X 10-25. 

Under the condition of “isentropic particle creation”, 
the changin g ra te of thermodynamic parameters is rep¬ 
resented by [13 


h 

n 


f 

T 


-( 0 -r), 


(26) 

(27) 


So for the selected P, the particle number density and 
the temperature of the cosmic fluid are given by 


n = Up 


T = Tp 



-3(l-/3) 


-3(1-/3)(7-1) 


(28) 

(29) 


would remain constant throughout the time evolution of 
the universe, i.e. the particle creation rate P tends to 
zero. Also we get radiation like solution for T oc R~^ as 
/3 —>■ 0 and 7 —>■ 4/3. 


4. CONCLUSION 


In this work, we studied the role of bulk viscosity in 
the framework of the particle creation mechanism in flat 
FRW universe and obtained analytical solutions through 
the full causal theory of Israel-Stewart. We focused on 
“isentropic particle production” process in which entropy 
per particle is constant. We solved the transport equa¬ 
tion of bulk viscosity for a linear form of T{H) adopted 
from CCDM scenario and obtained an expression for the 
bulk viscous pressure which is different from its coun¬ 
terpart equation in the first-order approach by an extra 
term that becomes important for very large R and may be 
responsible for late time acceleration. Using the derived 
bulk viscous pressure and a varying equation of state pa¬ 
rameter 7 ( 7 ?), we obtained Hubble parameter, particle 
number density and the temperature of the cosmic fluid 
for inflation and radiation eras. 

It is interesting that both p and H start with a finite 
value at t = 0 and end up with a small finite value at very 
large t. Therefore, in this picture the viscous model of¬ 
fers a desirable possibility for replacing the cosmological 
constant and explaining the late time acceleration of the 
universe. Further, It is notable that the inflationary solu¬ 
tion can be found with or without particle creation which 
is due to the constant bulk viscous coefficient. Also for 
later times, Hubble parameter has power-law behavior 
which corresponds to the radiation-dominated solutions. 
We showed that the solutions for the Hubble parameter 
of both epochs have good consistency with the equiva¬ 
lent solutions of ACDM model based on reported data 
in recent Planck collaboration papers. Using these data, 
we constrained the free parameter of particle creation in 
our model a < fldm = 0.308 ± 0.012. 

It is notable that N is an increasing function of time 
until a ^ 0.154, then it would remain constant forever. 
Also we found a radiation-like solution for the tempera¬ 
ture of the cosmic fluid as T oc R~^ for a ^ 0.154 and 
7 ^ 4/3. 
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